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Abstract. It is shown that the distribution functions of the diffusion coefficient are 
C^N | very similar in the standard model of quantum diffusion in a disordered metal and in 

a model of classical diffusion in a disordered medium: in both cases the distribution 
functions have lognormal tails, their part increasing with the increase of the disor- 
der. The similarity is based on a similar behaviour of the high-gradient operators 
pL^ ' determining the high-order cumulants. The one-loop renormalization-group correc- 

ly-j ■ tions make the anomalous dimension of the operator that governs the s-th cumulant 

proportional to s(s — 1) thus overtaking for large s the negative normal dimension. 
As behaviour of the ensemble-averaged diffusion coefficient is quite different in these 
models, it suggests that a possible universality in the distribution functions is inde- 
^ . pendent of the behaviour of average quantities. 

(N '■ 

: 1 Introduction 

ON . 

The absence of self-averaging in mesoscopic systems (see recent collections of review pa- 
pers [|1], |2j) put forward a question about probability distributions of different physical 
quantities. The conventional approach to systems with quenched disorder was based on 
describing measurements in a single sample by the ensemble-averaged quantities, i.e. the 
quantities averaged over all realizations of the quenched disorder. Discovery of the meso- 
scopic fluctuations (first of all, of the "universal" conductance fluctuations) manifested that 



this approach failed for the mesoscopic systems, at a length scale quite large as compared 
to a microscopic one. Instead, a statistical approach based on evaluating the distribution 
function of any physical quantity should be implemented. It could be reduced to evaluating 
the average value and the variance only, if a distribution is normal. However, a question 
about the shape of the distribution could not be solved a priory. 

For the problem of quantum diffusion, the mesoscopic scale is defined by either the ther- 
mal length or the phase-breaking length due to inelastic scatterings, both diverging when 
temperature goes to zero. Then, at T = the statistical approach is required for describ- 
ing any physical quantity at any scale. Such an approach, based on the renormalization- 
group (RG) analysis in the framework of the field-theoretical (cr model) description of the 

* Dedicated to Professor Hans A. Weidenmiiller on the occasion of his 60th birthday 
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quantum-diffusion problem , has been constructed for different physical quantities at 
the metallic side of the Anderson metal-insulator transition, and extended by means of 
the 2 + e expansion (near the critical dimensionality d — 2) for describing these quanti- 
ties in the vicinity of the transition (see ref. for reviews). A complete description of 
the transition itself should be in terms of changing the distribution of conductance rather 
than of a critical change of the average conductance only, as all irreducible moments (i.e. 
cumulants) of the conductance are of the same order in the region of the transition so 
that the conductance distribution is completely non-Gaussian. However, in this region the 
2 + e expansion provides quite a poor (at d — 3) description even for the average conduc- 
tance, as the four-loop calculations by Wegner|| proved. In view of this, is it possible that 
the statistical approach based on the 2 + e expansion could provide at least a qualitative 
description of changing the distribution function in approaching the transition? 

I will argue here that such a possibility does exist. A main qualitative prediction 
of the statistical approach is the appearance of the lognormal tails in the distribution 
functions of different physical quantities, their part increasing with the increase of the 
disorder. For some local quantities the distributions become completely lognormal still on 
the metal side of the transition |5[]. It provides a basis for the conjecture that statistically 
the metal-insulator transition reveals itself as a crossover from the normal to the lognormal 
distribution. Such a conjecture seems to be confirmed by recent numerical simulations || 
and by the analytical continuation of the one-loop RG results for different distribution 
functions from d — 2 + e to d — 1 that reproduces, after substituting a l-d value for 
the mean conductance, exact l-d results for these functions. It appears, therefore, that 
although the 2 + e expansion could not provide the reliable values for critical exponents, it 
could describe relations between the cumulants of the same physical quantities. The reason 
is that these relations could be independent of those particular properties of the model 
which determine behaviour of the ensemble-averaged quantities (e.g. the conductance). 
An important property of the nonlinear a model is the asymptotic freedom. The coupling 
constant (inversely proportional to the conductance) increases making the perturbative 
approach inapplicable in the region of the transition. The behaviour of the cumulants, 
however, is determined essentially by some additional operators. For each cumulant, the 
whole set of the additional operators is generated under the RG transformations. The 
eigenvalues of the RG equations that determine the critical exponents for the cumulants 
are governed by the operations of the group of permutations acting in this set. These 
operations seem independent of the renormalization of the coupling constant. 
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To provide arguments in favour of this viewpoint, I will consider a model of classical 
diffusion in a medium with the quenched disorder. The coupling constant of this model 
is not renormalized to all orders of the loop expansion. It does not make the model to 
be trivial, though, as the long-time diffusion is anomalous and characterized by the mean- 
square displacement that increases with time t slower than t. Although the behaviour of 
the coupling constant has nothing to do with the asymptotic freedom characteristic of the 
quantum-diffusion problem, the set of the eigenvalues of the RG equations governing the 
cumulants of the diffusion coefficients turns out to coincide with that for the quantum- 
diffusion problem. In both cases, the additional operators that govern the cumulants are 
irrelevant within the naive scaling but become relevant allowing for the one-loop RG cor- 
rections which give rise to nontrivial deviations of the distribution functions from Gaussian. 
It suggests that such a behaviour could be "universal" for very different models which by 
no means belong to the same universality class. As the classical-diffusion model is much 
simpler than the quantum one, a further analysis would be possible which is quite cumber- 
some in the quantum case. Due to this similarity, it could shed light on some properties 
of the Anderson transition. 

Before considering the classical-diffusion model, I will describe the principal results on 
the conductance distribution in a disordered metal. 

2 Conductance distribution in the quantum diffusion 
problem 

For the sake of comparison with the classical diffusion problem, I will outline a way of 
calculating the moments of conductance distribution using the field-theoretical technique. 
A detailed description can be found in ref.|J. 

The ensemble-averaged conductance g can be represented directly in terms of the gen- 
eralized nonlinear a model as follows: 



where L is size of the system, a and (5 are the vector indices in the <i-dimensional space, the 
brackets (. . .) stand for the averaging over all realizations of the disorder, and g is chosen 
to be dimensionless quantity measured in the units of e 2 /27r 2 h. The generating functional 
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is given by 

zlA] Iew{-F[Q;A]}VQ 
J exp{-F[Q}}VQ 

with F[Q] = F[Q; A = 0] being the functional of the standard nonlinear a modeljs], || 

F[Q} = ijd d rTi(dQ) 2 . (3) 

Here d stands for a usual gradient, and the coupling constant t = (7tzaD/8) _1 is related 
to g = (g xx ) as t = l6irL d ~ 2 /~g where D = v F l/d is the diffusion coefficient, v is the one- 
electron density of states, and I is the length of the mean free path. The 2N x 2N matrix 
Q obeys the constraints 

Q 2 = I, Tr Q = (4) 

and the integration in Eq.(^j) is carried out over independent elements of Q. In deriving the 
above field-theoretical representation from the initial model of the non-interacting electrons 
in the random potential, the replica trick has been used for the ensemble averaging so that 
the limit N = O should be taken on calculating the conductance (|l|). Alternatively, the 



super symmetric formulation [M, |TTJ could be used for calculating both the average value 



of conductance and its higher moments. 

The functional F[Q; A] in Eq.(0)is obtained from F[Q] by substituting the gradient dQ 
by the covariant derivative VQ: 

dQ^VQ = d a Q - [A a , Q] , (5) 

where the source field A is introduced to provide the direct representation of the conduc- 
tance in Eq.(|l]). 

A generalization of Eq.(|I|) for calculating the s-th cumulant of the conductance is 
straightforward |J: 



A=0 
iV=0 



(6) 



However, a further extension of the effective field-theoretical model was required in Eq.(^) 
to obtain all relevant contribution to the cumulants. So the generating functional Z[A] 
of Eq.(||) has been replaced by Z[A] obtained by substituting F[Q, A] for F[Q, A] into 
Eq.(Q) where the functional F[Q, A] of the extended nonlinear a model reads 

oo 

F[Q, A] = F[Q, A]+]TF S [Q, A] 

s=2 
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. 2n 

F a [Q, A] = z s S ai _ a2a d d rTr(T[ V ai Q) . (7) 

J i=i 

The bare values of the charges z s are inversely proportional to the bare value of the coupling 
constant t. The tensor 5 ai a2a in this expression arises from angular integration of the 
product n ai . . .n a% (n = p/p F ) so that it is proportional to the sum of all the possible 
products of the Kronecker symbols. 

The functional fl7|) alongside with the functional of the standard a model (Q) has been 
derived microscopically from the initial model of free electrons in the random potential 
V(r) chosen to be Gaussian with zero average and the short-range correlator 

(V(r)V(r')) = -J-*(r - r') 

2-KVT 

where r = l/v F is the mean time between elastic scatterings. 

The conductance cumulants (||) have been calculated with the renormalized functional 
Eqs.(HD, (0). It is important that the latter functional is not closed under the RG procedure, 
and may be renormalizable only in a larger space which includes for any integer s > 2 
(where 2s is the total number of VQ) all possible vertices of the following type: 

F g [ Q; A] = z* Jd d r [Tr(V Ql QV Q2 Q)] Sl . . . [Tr(V A Q . . . Vp 2m Q)] *\ . . (8) 

Here a set of integers s = (s\, . . . , s m , . . .) obeys the constraint 

ms m = s, (9) 

m>l 

and vertices with all possible partitions into the product of traces and all possible trans- 
positions of the vector indices and fa (each of them being repeated twice) are included. 
The notation (^j) is somewhat symbolic as it means that there are s m traces of the length 
2m while their vector- indices structures are different. 

The bare values of charges attached to all the additional vertices are equal zero. How- 
ever, they are not less "physical" then the original ones derived microscopically, Eq. ([?]). 
Thus, the vertices which made direct contribution to the cumulants @ are those with the 
set Si — s (and s m = for m/1) rather than the original ones (with the set s s — 1 and 
s m = for m 7^ s). 

All the high-gradient operators are naively irrelevant as the couplings z s — ► X d ~ 2s z s after 
the scaling XL — > L. However, the one-loop RG analysis which reduced, after deriving the 
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RG equations for the whole set of the couplings z s , to the analysis in terms of the group 
of permutations |5], |j gives the dimension a s of the couplings as follows 

a s = d-2s + es(s- 1) (10) 

where only the largest eigenvalue s(s — 1) of the RG transformations has been taken into 
account. At d = 2 one substitutes e by g~ x in this equation. 

There are two different types of contributions to the renormalized values of the con- 
ductance cumulants. The first one is obtained by the repeated differentiation in Eq.(^) of 
the generating functional Z, which includes only the standard functional of the nonlinear 
a model ([3]) with the substitution (|5|). The contributions of the second type are obtained 
with allowance for the high gradient vertices Eqs.([7j) and (|8|). On keeping the most impor- 
tant contribution of the second type (which is made for the s-th cumulant by the vertices 
with 2s gradient operators), one finds || omitting the vector indices that 

{ g 1 ' 8 , s^n (a) 

9 s ) ~ (11) 

Here go 3> 1 is the bare value of the conductance, 

u = \n— , (12) 
a 

o"o and a are the bare and renormalized value of the conductivity (not conductance). The 
number n Q that separates "high" and "low" moments turns out to be quite large in the 
region of validity of the RG approach, n ~ M _1 ln(L//), i.e. n ~ g Q in the region of the 
weak localization where g = go — ln(L/Z) ~ go ^> 1, and no ~ Qo/^iQo) on approaching 
the strong localization region where ~g drops down to 1. 

The exp(s 2 )-dependence of the moments is a generic feature of the lognormal distribu- 
tion as one can see from the following identity 



v s ex V (us 2 ) = (4?ra)1/2 J w s exp 



-^ln 2 ^ 

4m v . 



-■ (13) 
w 



As only the high moments of the distribution obey this dependence, only its tails are 
lognormal. In the weak-localization region, the low cumulants (|TT|a) are small as compared 
to variance (which is of order 1), and the bulk of the distribution is normal. On approaching 
the strong disorder region, the distribution becomes wide and non-Gaussian and the part 
of the lognormal tails increases. However, the small parameter (l/L) 2s in Eq.(|ll|b) which is 
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due to the normal part of the dimension fllPl ) of the couplings z s , prevents the tails becoming 
more essential at the metallic side of the transition. In calculating the cumulants of some 
other quantities, like a local density of states (DoS), the negative normal dimension of 
the high-gradient operators is exactly compensated by the positive normal dimension of 
the operators coupled to the cumulants of DoS so that the behaviour of the cumulants is 
governed completely by the anomalous dimension. In this case, the distribution function 
becomes completely lognormal when u ~ 1 which is still in the metal regime. || 

The distribution of the global DoS is similar to that of the conductance. Thus, due to 
the Einstein relation, the distribution of the diffusion coefficient (which, for the quantum 
diffusion problem, has not been directly calculated) has similar features. I will show below 
that the same behaviour of the distribution of the diffusion coefficient is characteristic of 
some model for classical diffusion in disordered media. 



3 Classical diffusion in media with the quenched dis- 
order 

Diffusion of classical particles in random media has attracted considerable attention (see 
ref . [Il2| for a recent review) . A set of mo dels@-0) is given by the Langevin equation that 
describes random walks of a particle under the influence of a thermal noise rj(t) in the 
medium with the quenched disorder governed by a random-drift field v(r): 



r = v (r) + rj(t) 

The thermal noise r\ is chosen to be Gaussian white noise with a zero average and 



(14) 



V*{t)vp(t') = 2DS aP S(t-t'), 



(15) 



where D is the diffusion coefficient. 

Equivalently, averaging over the random noise one describes these random walks in 
terms of the Fokker-Planck (FP) equation 



d_ 

di 



+ d a (v a - Dd a 



p(r,t) = 



(16) 



with p(r, t) being the local density of the diffusing particles. 



The exact solution obtained by Sinai |1| to the model Eqs. (|T4])-(|T6|) in one dimension 
(with the Gaussian short-range correlator for the random-drift fields v ) has revealed the 
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anomalous diffusion manifested by the strongly sub-diffusive long-time behaviour of the 
random walks. Later such a diffusion has been considered as a possible universal source 



for the excess (1//) current noise [18 



However, it appeared that anomalous diffusion occurs only for d < 2, at least in the 
weak-disorder d = 2 turned out to be the upper critical dimensionality, and the 

long-time behaviour of the random walks could be super-diffusive as well as sub- diffusive, 
depending on some features of the model. The point is that a vector character of the 
quenched disorder in d > 1 gives more freedom in choosing correlator of the fields v. 
In general, assuming the Gaussian nature of the random-drift fields, i.e. irrelevancy of 
the higher-order correlations (the assumption that could be proved in a subsequent RG 
analysis), and choosing (v) = 0, one may write has 

(v a {r)v p {r')) = lo F af} (r-r') (17) 

where 70 characterizes the strength of the disorder. 



The simplest possible choice of the correlator ]14 

Fa/3 = S aP 6(r - r) (18) 

leads to the model where corrections to the effective diffusion coefficient D(t) vanish in 
the long-time limit. Two other specific models for F a p arise naturally£l _ @ ) on imposing 
constraints on the random field v choosing it either potential (v a = <9 a $) or solenoidal 
(d a v a (r) = 0). For these models the Fourier transform of the correlator ([17]) reads 



a/3 



a J 3 , potential field (a) 

(19) 

$a0 — ^—J- , solenoidal field (b) 



The RG analysis has shownEl^B) that at d — 2 the potential disorder leads to the sub- 
diffusive behaviour in the long-time limit while the solenoidal disorder leads to the super- 
diffusive one: 

a , v f D (t/r)- g °/ 2 , potential field (a) 
D(t) = ^ (r 2 (t)\ oc <^ (20) 
dt X 1 { D \n 1/2 (t/r), solenoidal field (6) 

Here 

— ^ (2D 



4nD 
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is the bare value of the coupling constant g (dimensionless at d — 2) which plays a role 
of the effective weak-disorder parameter, Do is the bare value of the diffusion coefficient 
(describing diffusion at t ~ r), and r is some small time which provides the ultraviolet 
cutoff necessary to regularize the model. 

Note that the results of Eq.(^) are asymptotically exact in the weak-disorder limit. In 
the case of the solenoidal disorder, the coupling g is decreased by the RG transformation 
from its bare value go 1 (the "zero-charge" situation) which makes the strong-disorder 
region unattainable. In the case of the potential disorder, there is no renormalization of 
the coupling constant g at all. It has been firstly demonstrated up to the two-loop order 
of the RG analysis but has later been proved [[OJ to be perturbatively exact: the 



renormalization is absent in all the orders of the loop expansion. It does not make the 
model trivial as both the strength of the disorder 7 and the diffusion coefficient D (as well 
as other physical quantities, such as mobility) are renormalized which leads, in particular, 
to the sub- diffusion, Eq.(|20|a). 

The quenched disorder may enter the FP equation (|T6|) also via the spatial variation 
of the diffusion coefficient -D(r) = D + 25D(r) where D is the uniform part of D{j). One 
can take the limit v = and define the disorder in terms of the distribution P(D) of the 
spatially uncorrected diffusion coefficients. The distribution is governed by its cumulants 



<fD(n)<fD(r 2 ))) = 2!^( ri -r 2 ), 

(22) 

5D(y l )...5D(r s )\ = s\ g {s) 5{v x - r 2 ) . . . <f(n - r a ), 



where the additional couplings could have changed under the renormalization. 

Such a model was believed to be trivial in the weak-disorder limit as the operators 
described the quenched fluctuations of the diffusion coefficient are irrelevant under the RG 
transformations W% at any dimensionality d > so that the long-time diffusive behaviour 
remains normal and the limiting distribution of the diffusion coefficients turns out to be 
Gaussian whatever is the initial choice of the constants g^ B >. However, as it has been 
recently suggested by Tsai and Shapir [P0| , this model could be analyzed in terms of the 
d = + e expansion which could shed some light at the diffusive behaviour in the strong- 
disorder limit. In the case of the annealed disorder, the analysis of ref. |21J has shown that 
in the "strong-disorder" limit (i.e. for some nonzero e) a nontrivial fixed distribution of the 
diffusion coefficient existed which turned out to be the inverse Gaussian one. Unfortunately, 
it appeared to be impossible |2(J to have extended the results of this analysis to the more 
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interesting case of the quenched disorder that corresponds to the N = replica limit. (The 
annealed limit corresponds to the N = 1 case, i.e. to the averaging of the non-replicated 
functional). 

Notwithstanding a viability of + e expansion for describing the strong disorder in 
any physical dimensionality and the lack of analysis in the case of the quenched disorder, 
the appearance of the nontrivial fixed point |2(J seems to be quite important. However, 
the "zero-dimensional" model defined by the FP equation (|T6|) with v = and by the 
distribution ( p2[ ) has some quite specific features. The most important is that all the 
operators coupled to the diffusion cumulants become marginal at the critical dimensionality 
d = (i.e. the naive scaling dimension for each operator becomes equal zero). In any model 
with d 0, the higher the order of cumulant, the more irrelevant becomes the appropriate 
operator. Thus, in the model of the quantum diffusion described in the Section @ the 
naive dimension of the 2s-th high-gradient operator (which defines the s-th cumulant of the 
conductance and, to some extent, of the diffusion coefficient) decreases proportionally to 
— s as in Eq . (|10|) . It makes impossible an impact of the operators with the higher number 
of gradients on those with the lower one that, in turn, makes implausible the appearance 
of a new fixed point for the distribution (at least, starting from the weak-disorder limit 
where the dimensional analysis has some sense). It does not mean, however, that the 
distribution remains trivial (i.e. normal). As was shown in the Section (0), to the one-loop 
order the dimension of the high-gradient operators becomes proportional to s 2 so that they 
could become relevant. It leads to the increasing ( |TTD of the higher cumulants which define 
nontrivial tails of the distribution. As these tails become more essential on approaching the 
strong disorder regime, the shape of the distribution could have changed completely in the 
strong-disorder limit || . As all the cumulants have different critical exponents (which are 
not proportional to the order of the cumulant), the change in the shape of the distribution 
has a character of a smooth crossover that is not associated with a new fixed point. 

The goal of further considerations is to show that a similar scenario can take place in 
a classical-diffusion model. I will consider the FP equation ( |16|) that includes both the 
random drifts v(r) defined with the correlator Eqs. (p~7|)— (|19|) , and the random diffusion 
coefficients defined with the cumulants fl22l) . 

Before proceeding with the renormalization, I will describe a possible realization of the 
random-drifts models. 
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3.1 RELATION TO LATTICE HOPPING MODELS 

The lattice hopping model which goes over to the FP equation ( ]T6|) in the continuum limit 
is defined by the master equation for the probability p r {t) of finding the particle at the 
r-th site at the moment t: 

^ = Y,( W rr'Pr'-W r , rPr ). (23) 
r' 

Here W rr r is the probability of hopping from site r' to site r which slightly fluctuates around 
the hopping probability W° in a regular lattice: 

W rr , = W?_ r , + 5W rr , } (24) 

with 5W rr > describing the quenched weak disorder on the lattice. 

The parameters of the continuum model described by the FP equation (|T6| ) are readily 
expressed in terms of the hopping probabilities: 

D(r) = 1 J> - v'fW TT > (25) 

r' 

v(r) = 5>-r')W rr , (26) 



The random drift velocity given by Eq. fl26|) vanishes in all realizations of the disorder 
until the asymmetry in the hopping {8W TT > 7^ 5W r > r ) is included. Therefore, the symmetric 
hopping model goes over to the FP equations with v = where the disorder in D described 
by the cumulants (|22"D reflects the symmetric disorder in W. In this case, the analysis of 
the continuum limit showed the weak disorder to be irrelevant for any d. 

The asymmetry could be due to the presence of magnetic or charged impurities which 
leads to the continuum model with either the solenoidal or the potential disorder. Following 
ref . [ 1 7 1 , 1 will show that in the presence of a random electric field E(r) of charged impurities 



the continuum- limit model includes the potential random drifts (|I~9|a). If E is the only 
source of randomness, the probability of thermally activated hops at a distance b is 



Wr, T+b = exp 



eE(r)b 

2kT 



(27) 



Assuming then a random distribution of the charged impurities on the lattice and a global 
electro-neutrality, and representing the Fourier transform of E(r) as 

E(q) = }^ e i ex P(«q • r j)> 
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with ej = ±e being the charges of impurities with random coordinates rj and e being the 
dielectric constant, one finds 

(EM)EM)} = <^M! V ( 28 ) 

£ q 

where C is the concentration of the charged impurities. In the continuum limit, using 
the linear in E approximation (5W oc E) and assuming for simplicity the nearest-neighbor 
hopping only, one arrives at the potential-random-drifts model with the correlator Eqs.flHj), 
(|HS|). The strength of disorder 7 in Eq.(|I7|) is related to the parameters of the hopping 
model as 

ttVC [ W za 2 \ 2 ttVC 2 

10 ~^{kfy\^~ ) = ^ry D ° (29) 

where z is the coordination number of the lattice, a is the lattice constant, Wo = is 
a regular part of the nearest-neighbor hopping probability, and Dq is the bare value of 
the diffusion coefficient. The randomness in W results also in some distribution P{D) of 
the diffusion coefficient whose cumulants (p2^ ) could be directly expressed in terms of the 
cumulants of 5W. 

In real lattices the correlator (p8|) becomes non-singular (oc q a qp) for q ^> r where r is 
a screening radius. The continuum model with the correlator (|19|a) could then describe only 
the random walks at a distance not exceeding r . Some models where r$ is macroscopically 



large are described in ref.|T7|]. In general, r takes the part of the infrared cutoff similar to 
that of the phase-breaking length in the weak-localization theory and defines the mesoscopic 
scale for the classical-diffusion problem. Below I will discuss the continuum model with 
the unscreened correlator (|I§D. 

3.2 THE EFFECTIVE FUNCTIONAL 

Random walks of the particle started at the initial moment t — at the point r is charac- 
terized by the probability to find it at the point r' in time t 

G(r',r;t) = (p(r',t)p(r,0)) (30) 

which is the Green's function of the FP equation (|T^|). The asymptotic properties are 
governed by the long-time behaviour of this function, or by low-frequency behaviour of 
its Fourier transform, G(r,r';cu). To perform the ensemble averaging, it is convenient to 
represent the latter as a functional integral over the conjugate complex fields y?(r) and 
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G(r ' r ' w) " JeWriTnpVtp ' (31) 



where the effective action functional is given by 

S\(p, (f] = Jd d r[i7puj(p + v a (d(Jp)<p — Dd a lpd a (p\ (32) 

The averaging is performed by the standard replica trick: the fields </?(r) and ip(r) and 
the functional integration in Eq. (|3~I]) are iV-replicated and the independent averaging over 
the numerator and denominator in Eq.(|3lD is justified in the replica limit N = (that 
should be taken in the final results). Taking into account the disorder both in the random 
drifts v and in the diffusion coefficients D which is defined by Eqs. (p]) and (|22|), one 
deduces the effective action 

S[jp~, ip] = {S + S int + S cum } [Tp, tp\ (33) 
which should be substituted for that given by Eq. (|32"D into Eq. (|3l"l) where the functional 



integration should be performed over all components of the fields ip = (fx, . . . ,(fN and 
lp = ip 1 ,...,ip N . Here 



S [lp,<p] = Jd d rlp(iuj + Dd 2 )<p (34) 
S^&M = y Jd d rd d r\d a Tpv) r F a p{v-r>){dplpip) rl (35) 

oo . s 

Scum&,¥>} = iJ29 (s) S(sy, S {s) = d d r ]J{d ai Tpd ai (p) r (36) 

^ — O <J A 1 



=2 " i=l 



The RG analysis of the functional (34]), (^5|) in the upper critical dimensionality d = 2 



has been used for describing the anomalous long-time behaviour of the average diffusion 
coefficient, Eq.(pO|). 

The higher-order gradient operators, Eq . (|36|) , describe the renormalization of the cu- 
mulants ( p2|) of the distribution P(D). Under rescaling L — > XL, — > \^ s ^ d g^ so that 
the naive dimension of behaves like that for the high-gradient operators in the quantum 
diffusion problem, and are irrelevant. I will show now that in the case of the potential 
random drifts, Eq. (|l9|a) , the one-loop RG corrections overturn this conclusion and make 
the scaling dimensions of the high-order gradient operators positive. 
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3.3 ONE-LOOP RENORMALIZATION 



The renormalization is performed by expanding exp(iS int + iS cum ) in a power series, in- 
tegrating with the weight exp(2iS>o[^ ' Vol) over the "fast" components of the fields (this 
integration will be denoted as (. . .) ), and exponentiating the results of the integration. 
Here (p(r) (and <^(r)) is decomposed into the sum of "slow", <£>(r), and "fast", ip (r) 
components where 

£(r) = £ ^(q)e^ r , ^ (r) = £ cp(q)e^ r , 

q<Mo Mo<q<Qo 

< A < 1 is the scaling parameter, qo is the ultraviolet cutoff, e.g. the inverse lattice 
constant in the lattice realization of the model. In the one-loop order (i.e. in the first order 
in powers of the coupling constant g, Eq. (|2lD ), the logarithmic (at d = 2) contribution to 
g(s) j s made only by the term 

- (S int S (s) ) o = -gAR{S {s) }. (37) 

Here the one-loop RG operator R is introduced, A = In A -1 . The relevant contractions 
(i.e. those which give the logarithmic term) have the structure 

(^^o) (rtdf&g ) Q , (38) 

where fields (p are taken from S int , and fields Tp Q from St s y As (^o</? ) is j us t a diffusion 
propagator, the contractions containing more gradients are irrelevant. That is why there 
is no contributions from the terms like (S^ um ) Q ,and no feedback from the terms with the 
higher number of gradients to those with the lower one. 

It is evident from Eq.(^) that the index structure of the action (|36|) is not conserved 
under the renormalization. A set of additional operators is generated having the structure 

d a Tp a d^ b d^ c d 5 v d ... (39) 

with all possible permutations of the vector indices a, (3, . . . and the replica indices a,b, . . ., 
each index being repeated twice, where the summation over repeated replica indices from 

1 to N and over repeated vector indices from 1 to d is implied. These operators are 
"unphysical" in a sense that the distribution P(D) is defined by the renormalization of the 
initial operators ([36]) only. Then, only those operators are of interest which have the RG 
feedback to the initial ones. 
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To classify all the additional operators one introduces matrix notations 

Q ab = d a Tp a d a ip\ V ab = d a <p a d a <p\ V ab = d a Tp a d a Tp b . (40) 



In these notations, the contractions Q38|) could involve only the matrices Q and V. On 
calculating the action of the RG operator ( P7|) on the product QQ one finds 

RlQ^Q^X = \A(Q ab Q cd + Q ad Q cb ) + (2 - \A\)V ac V bd \ (41) 

Here A discriminates the models described above: A = 1 for the potential random drifts, 
Eq. (|19^t) , A = — 1 for the solenoidal ones, Eq.(|T9|b), and A = for the spatially uncorre- 
latted random drifts, Eq. (|i~8|) . 

As R does not act on the matrices V (they do not contain the field Tp to be contracted 
with the field (p from the Si nt ), their number in any operator could not be reduced under the 
renormalization. Therefore, the RG equations have a triangular structure: the operators 
containing only the matrices Q are not influenced on by those containing VV (the number 
of V is evidently equal to the number of V in any operator). 

The initial operator (^) is proportional to [Tr Q] s , so that one can restrict the RG 
analysis to the subset of the operators containing only the matrices Q. There is no one- 
loop corrections to in the case of the uncorrelated random drifts, Eq.(|l8|), as only 
the VV term arises on renormalizing the initial action (p6|). I consider the case of the 
potential random drifts, Eq.(|l9"la). The action (f|l|) of the RG operator R on the chosen 
subset reduces to simple permutations of the matrix indices of Q. It can be analyzed in 
terms of the group of permutations, similar to the case of quantum diffusion [^, ||. Any 
operator containing s matrices Q (i.e. 2s gradients of the fields Tp and </?) may be written 
down in a form similar to that for 2s-gradient operators in the quantum-diffusion problem 
(compare to Eq.(|)): 

Sf[Q]=g® Jd d r {[Tt(Q)] Si ...[Tr(Q) m p . . .} (42) 

Here the set of integers s obeys the constraint ([]), the initial operator (|36D corresponding 
to s = (0, . . . , 0, s). The renormalization of the s-th cumulants results from solving the RG 
equations for the whole set of g^, the bare values of all the additional charges being equal 
to zero. 

The RG equations may be written down as follows 

fU«l?> (43) 
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where is represented as some ket- vector defined by the "occupation numbers" s m 

g® = \s) = \1 S1 2 S2 . . .m s ™ . . .), (44) 
and R is expressed in terms of the "creation" and "annihilation" Bose operators at, a 

Z ^ m>l m>l 

+ J2 [i m + a \n a \ a m+i + ml a} m+l a m ai L (45) 
m,l>l ' 

these operators acting on the states \s) as follows: 

a m |l Sl . . .m Sm . . .) = s m \l Sl . . .m Sm ~ l . . .) 

4.1 v (46) 
a^|l Sl • • • m Sm . . .) = |l Sl . . . m Sm+L . . .} 

So, the RG transformation is expressed in the form of some one-dimensional field theory 
with the "hamiltonian" (fl~5|). The first line in Eq . (|45|) represents the trivial diagonal part in 
the operation (f4l|) while the second one represents the transposition between indices of Q in 
Eq. (fUD . The hamiltonian (f4"5|) differs only by the trivial diagonal part from that derived for 
representing the RG transformations for the composite operators of the unitary nonlinear 
a model that describes the quantum diffusion problem in the presence of a magnetic field. 

The operator (|45|) which is just a representation of the simple symmetric group opera- 
tion (all possible pair transpositions ( fH"D in a set of 2s pairs of indices) can be diagonalized 
exactly §]. The eigenvectors are given by 

\p) = J29(s)x p (s)\s) (47) 
{*} 

where the summation is performed over all the partitions s = si, . . . , s m . . . of s obeying 
the constraint @, g(s) = s\/l\ m m Sm s m \ is the number of elements in the class defined by 
the partition s, and Xp{s) are the characters of irreducible representation of the symmetric 
group characterized by the Young frame p having boxes of length px, . . . , p rn where J2 m Pm = 
s. The appropriate eigenvalues are given by 

a , (p -) = f<i^ + £ ^- 2 '" + 1 > (48) 



To verify Eqs. (|47|) and fl48|) , it is necessary to remember that the operator ( ^q) is non- 
Hermitian, and and a defined by Eq. (|46|) are not mutually conjugate. In view of this, 
the set of bra-vectors orthogonal to the ket-vectors (|47|) is given by 



(Pi = (1/s\) J2x p (sM. (49) 



16 



The maximum eigenvalue corresponds to the eigenvector characterized by the one-line 
Young frame with p\ = s, p m = for m > 1 for which xp{$) = 1 f° r an so that it is equal 
to s(s — 1), as in the case of the quantum diffusion. Note that it could be verified without 
any reference to the representations of the permutation group, acting by the transverse 
operator R on the bra- vector ( fi9|) which is proportional in this case just to XX s|- Thus, 
with the one-loop accuracy, the dimension of the operators coupled to the moments of the 
diffusion coefficient is given by 

a s = -(s-l)d + g s(s-l), (50) 

very similar to the case of the quantum diffusion, Eq. (|T0|) . 

Note that the minimum eigenvalue in Eq. (fI8D (which corresponds to the eigenvector 
characterized by the one-column Young frame) equals zero. As results from Eq. ([4lD , the 
eigenvalues for the solenoidal random-drifts differ in sign from those for the potential case. 
Thus, there is no positive exponents in the solenoidal case so that the dimensions of the 
operators coupled to the diffusion cumulants remain negative, and the distribution scales 
to the normal one. Therefore, of the three random-drifts model considered, the distribution 
becomes nontrivial and similar to the quantum-diffusion case only for the potential drifts. 



4 Comparison of the results for the quantum and 
classical diffusion 

There are two types of contributions into the fluctuations of the diffusion coefficient in 
the random- walks model considered, similar to the quantum diffusion problem, Eq. (|TTD . 
The "normal" one is given only by the functional (p5|). It diverges in the infrared limit 
thus making the fractional fluctuation (((5D) 2 )) / D 2 oc gQ to be independent of the size 
of the system, analogous to the universal conductance fluctuation (the constant go has 
nonuniversal dependence on the disorder parameter, Eq . (|2l|) , which restricts the analogy). 
This contribution has been considered in detail in ref. fl7| and is of no interest for the present 



considerations. The additional contribution to the diffusion cumulants ( p2[) is governed by 
the dimensions of the couplings g s in the "additional" functional (|3(J). Keeping only the 
maximum eigenvalue, as in Eq.([5(]), one finds in the critical dimensionality d = 2 

/ j \ 2(s-l)-g s(s-l) 

(^r)oc m (si) 
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It is quite similar to the "additional" contribution to the conductance cumulants, Eq.(|TT|b), 
in the quantum-diffusion problem (as the appropriate contribution to the cumulants of the 
density of states has the same form ||, one can write down the same expressions directly 
for the quantum- diffusion cumulants). To make the analogy more striking, one substitutes 
into Eq.(|TT|p) the value of the parameter u, Eq. (|T2"D , in the weak-localization limit at d = 2, 
u = g^ 1 ln(L/Z) = u, which gives 

J'» <* (j) - (52) 

Note that both expressions (|5T|) and (|52|) are applicable only for the high-order cumulants 
as all the eigenvalues of the RG equations but the maximum ones have been neglected in 
the derivation. Nevertheless, even for small s > 2 these expressions give a qualitatively 
valid result that the low-order cumulants are small as compared to the variance so that 
the bulk of the distribution is Gaussian. The high moments (with s^g^ 1 for Eq.(|5"TD and 
s^g in Eq.fl5"2"D) govern, in accordance with the identity ([TBI), the lognormal tails of the 
distributions which have the following shape 



a 2 



f(SX) oc ^exp 



-—In 2 (6Xa 2 e- u 
Au v 



(53) 



An appropriate derivation is described in detail in ref . || . Here 5X stands for either SD or 
for g in the classical or quantum problems, respectively, a = L/l, and u is given above for 
the quantum problem, and equals goln(L/Z) for the classical one. 

In the classical case, u has only a trivial dependence on the scale given by the loga- 
rithm, as there is no renormalization of the coupling constant. In the quantum case, u 
gives just the weak-disorder limit for the parameter u = ln(<To/cx), Eq. (|12D . So the whole 
dependence on the renormalization of the coupling constant is absorbed by the parameter 
u which is substituted to the lognormal distribution, its appearance being governed by the 
renormalization of the high-gradient operators. That is why one can separate the renor- 
malization of the average values from that of the higher order cumulants. In such a way, 
one reproduces exact one-dimensional results for lognormal distributions by substituting 
the exact one-dimensional value of u into the formulae similar to (|53"1). It provides a basis 
for the conjecture that one can obtained a reasonable description of the distributions near 
the transition just by substituting exp(w) ~ \g — g c \ u with a proper choice of the critical 
exponent v. 
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5 Discussion 



It has been demonstrated that the shape of the distribution of the diffusion coefficient in 
the model of classical diffusion in a medium with quenched random drifts proved to be 
very similar for that of the conductance distribution in a weakly disordered metal. In both 
cases the distributions turn out to be almost Gaussian in the weak-disorder limit but have 
slowly decreasing lognormal tails, and the part of the tails increase with increasing the 
disorder. This similarity occurs although the coupling constants in the field-theoretical 
models describing the quantum and the classical diffusion behave very differently. The 
asymptotic freedom of the nonlinear a model, i.e. the increase of the coupling constant 
(inversely proportional to the conductance) with increasing a scale, is believed to govern 
the Anderson transition. No transition occurs in the classical diffusion problem described 
by the field theory with the (<900) 2 -type interaction where in the case of the potential 
disorder a perturbative renormalization of the coupling constant proves to be absent thus 
leading to the sub-diffusion, Eq . (|20|a) . The reason for the similarity is that the RG equa- 
tions governing the cumulants of the distributions in both cases are classified according to 
the same irreducible representations of the group of permutations. But deriving the RG 
equations for the high-gradient operators proved to be much easier in the classical-diffusion 
model. Then the higher-loop RG analysis of the high-gradient operators which could check 
a viability of the one-loop results seems hardly possible in the quantum model but could 
appear quite straightforward for the classical one. Thus studying the classical-diffusion 
problem described here gives a possibility to learn more about the quantum diffusion in 
disordered media. 
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